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Abstract 



X 

5_i ■ Several theorems are demonstrated that determine the sufficient 



conditions for the existence of synchronized states (periodical and 
chaotic) and also of travelling waves in a CML. Also are analytically 
proven the existence of period-doubling cascades for the mentioned 
patterns. The temporal state of any oscillators are completely charac- 
terized. The given results are valid for a number of arbitrary oscillators 
whose individual dynamics is ruled by an arbitrary C 2 function. 
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Systems showing patterns as a consequence of the interaction among their 
diverse components are really frequent, in any field that one can imagine: 
neuronal activity within the brain, or the function of organs as a whole 
within the body, drivers on a motorway, birds flying in a group, a network 
of computers, coupled lasers, crystal growth, etc. 

The result of the interaction of the individual elements generates struc- 
tures that manifest in the system as a whole. In these processes, one should 
consider two things: the behavior of any individual and the interaction among 
them. If we consider the traffic example, it is clear that the behavior of an 
individual driver, that is his decision to drive in a particular way or another, 
is certainly different when there are few cars on a motorway or when there 
is a traffic jam (in which case he will be guided by traffic patterns). 

Broadly speaking, all of these systems consist of a group of elements 
coupled by some kind of process, and at the same time, every element of the 
group is ruled by its own local dynamics. The understanding of such systems 
is extraordinarily complicated, since there are no particular mathematical 
tools developed to study them. One way to confront this problem is to 
discretize spatial and temporal variables as well as to fix inter-individual 
interactions as well as the individual dynamics. The result is a Coupled 
Map Lattice (CML) [T]: a chain of coupled elements (called oscillators), each 
situated on a discrete point of the lattice, whose individual dynamics is ruled 
by a discrete map. Despite the spatial and temporal variables are discretized, 
state variables remain continuous. 

In the last few years, CML have been extensively studied since the work 
of Kaneko and colaborators [21 EJ 0J El [6] , and from the beginning, they have 
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shown themselves to be exceptional modelling spatially extended systems. 
The use of this study has been extended into diverse scientific branches with 
an extraordinary variety of applications in physics, biology, chemistry, social 
sciences, and engineering modeling. [8] 

A typical evolution equation for a CML [TJ is given by 

m 

Xl ( n + 1) = (1 - a)f(Xi(n)) + - £ f(*j W) (1) 

3=1 

% = 1, m 

where Xi{n) represents the state of the oscilator located at node "i" of the 
lattice, in the instant "n". The parameter a weights the coupling among 
oscilators. Periodic conditions are assumed in the boundaries, given as 

Xi(n) = X i+m (n) Vz 

Depending on the value of a , the system behavior changes from the inde- 
pendent evolution of each oscilator (for a = 0) up to a mean field approach 
(for a = 1). For intermediate values < a < 1 the system is ruled by both 
local and global mechanisms. 

The general form of the coupling term is given by 

m 

m 3=1 
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where the measure the weights between the j-th oscilator and the i-th 
one. To achieve a symmetrical and spatially invariant coupling, it is usually 
taken Wij = w\i-j\. Sometimes, the coupling term will be written as 



(nearest-neighbor coupling). However, this last description is not adequate 
when we are dealing with a supercritical bifurcation threshold, because the 
coherence lengths are usually quite large [9]. Given that, in this paper, we 
want to study bifurcations in CML, we will use the mean field approach. 

Another important point, that must be considered, is the updating of os- 
cillators; they can be synchronous (all oscillators are updated simultaneously) 
or asynchronous (oscillators are updated one at a time) [lOl ITT] . Choosing 
one or the other depends on whether oscillators communicate among them 
much quicker that the updating time of the system as a whole, which is ruled 
by the evolution equation (Op). In this paper we will refer to synchronous 
systems. 

In the scientific literature, the majority of the results, referring to CML, 
are numerical results, as we will see later. The awesome richness of numer- 
ical results is restricted by a fixed and finite set of parameter values, and 
a finite number of oscillators in the CML, which supposes a limitation for 
adequate understanding of certain phenomena. In particular, the transition 



ni 




(mean field), or 



^[/(A i _ 1 (n)) + /(A i+1 (n))] 
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to chaos by period duplication needs the period to tend to infinity. It is also 
necessary the number of oscillators to be infinite, in a finite region, for the 
understanding of the onset of turbulence in fluids and plasmas; otherwise, 
there would be a cutoff in the wave numbers that could be studied because 
the lattice would have a finite spatial resolution. Mathematical proofs would 
be desirable to characterize syncronized states, traveller wave bifurcations 
and other behaviours. Fortunately, numerical results point out us what to 
look for and where. 

In this paper, analytical proofs, in CML, of the existence of syncronized 
states and travelling waves will be given. It will be proved that both pat- 
terns will go under a period doubling cascade as / , in ([TJ), does. These 
behaviours will be completely characterized, giving analytical expressions of 
the temporal evolution of every oscillator. 

The fixed points of CML, generated in period-doubling cascades, will be 
essentially the fixed points of f m2k (m number of oscillators in CML). As / 
determines the individual dynamics (see jTJ)), what is shown is the emergence 
of global properties from the local dynamics of a single oscilator. 

We have tried to keep the widest generality in the results; therefore, 
theorems have been proved using an arbtirary C 2 function f(x; r), instead of 
working with the logistic equation (or any toplogically conjugated functions) 
as usual. 

Perturbative methods will be used to obtain analytical solutions. The 
inversion of functional matrices of arbitrary size is fundamental in the proofs 
of the theorems; given that whenever the inverse matrix exists, it is unique, 
it will not be necessary to explain the calculation leading to it: it will be 
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enough to check that the proposed matrix (in the corresponding theorem) 
is the inverse matrix one was looking for. The matrices appearing during 
the demostration process will not be circulant; therefore, usual analytical 
inversion processes of circulant matrix inversion will not be valid. 

This paper is organized as follows. First, synchronized states will be 
considered, this solution being quite straightforward, it will indicate how to 
face up to the more complicated travelling waves in the next section. Both 
results will be used to study the period-doubling cascades of the patterns. 
The paper concludes with a section indicating connections of this work with 
other researchs. 

1 Regular and chaotic synchronization 

In this section straightforward analytical results will be presented for syn- 
chronization in CML, that is, for all the oscilators having the same value at 
anytime. This is a striking behaviour, in particular when chaotic syncroniza- 
tion is produced, where chaotic systems are very sensitive to perturbations 
and it is supposed that any slight modification generated by the coupling of 
the oscillators of CML would destroy the synchronization. The mathematical 
approach to this problem is far from being unique [T2] . 

Let 

m 

X t (n + 1) = (1 - a)f(X t (n)) + - V f(X t (n)) i = l,...,m (2) 

m. — ' 
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be the CML, with m oscillators, being a the coupling parameter and f(x) 
a function depending on a parameter r, in function of which the system 
Un+i = f(y n ] r ) shows fixed points for some arbitrary period p. 

1.1 Fixed points of the system. Stationary synchro- 
nized state 

It is straightforward to get the fixed points of the system. If the function 
f(x) has a fixed point in x* then (x*,x*, . . . m > ,x*) will be a fixed point of 
the system given by (J2J), since if 



Xi(n) = x* i = 1, . . . , m 



it turns out that 



f(X i {n)) = f(x*) = x* = X i {n) i = 1, . . . , 



m 



and 



rn 



X t {n + 1) = (1 - a)f(x*) + - J2f(x*) = f(x*) 



= X* 



i=l 



I = 



1,... 



m 



as it was wanted to prove. 



It is then deduced that: 



X{n) = (x* , x 



.* rn) * 



(3) 



, . . . ? 



is a stationary sychronized state of the system. 
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Figure 1: Stationary synchronized state. CML with f(x;r) = rx(l — x), 
r = 1.0, a = 0.1 and e = 0.1 

It is observed that if chosen x* andr for f(x*; r) to determine a periodical 
or chaotic evolution of x*, then the result would be that CML would have 
correspondingly periodical or chaotic synchronization, being this a proof of 
the existence of synchronized states both periodical and chaotic. In contrast 
from the CML with nearest neighbour coupling, there is not an upper limit 
in the number of oscillators "m" such that stable synchronous chaotic state 
exists [13]. See figures HHH 

Let us now study the linear stability of fixed points, where the eigenvalues 
of jacobian matrix will be calculated. 
The jacobian matrix is given by: 
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Figure 2: Period-2 synchronized state. CML with f(x;r) = rx(l — x), r 
3.2, a = 0.1 and e = 0.1 




Figure 3: Period-4 synchronized state. CML with f(x;r) = rx(l — x), r 
3.4985, a = 0.1 and e = 0.1 
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In 



dXjjn + 1) 
dXj(n) 



x* 



m—l 



a) 



a 
in 



\ 



a 
rn 



f'(x* 



where X* = ( *) and its eigenvalues are 



A = f(x*) single 
A = (1 — a)f'(x*) multiplicity (m — 1) 

therefore, the fixed point given by ((3l), or what would be the same, the 
stationary schronized state, would be stable whenever x* is a stable fixed 
point of f(x). 
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1.2 Periodical synchronized states 

The existence of periodical syncronized states is reflected in theorem 1, shown 
below. The way to proceed with this proof is similar to the one used to obtain 
the stationary synchronized state. 

Theorem 1. Let {xj, x%, x*} be ap-periodic orbit of C 1 function /. Then 
the CML given by 

m 

X^n + 1) = (1 - a)f(Xi(n)) + - V f{X t {n)) i = 1, . . . , m 

m ^-^ 

t=i 

i) shows a synchronized state of the same period as the function /. 
The sychronized states are as follows 

(r* r* m) r*V i 

ii) the sychronized states (x*,x*, . . , x*)j=i y ... iP have the same sta- 
bility as the fixed point x* of f p , j = 1, . . . , p. 

Proof 

i) Taking, for a given time n, 

Xi(n) — x\ i — 1, . . . , m 

results in 

F(Xi(n)) = Xi(n) i = l, ...,m 
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with the first iteration of CML being: 

m 

X t (n + 1) = (1 - a)f(xl) + fi<) = f«) < = 1, • • • , m 

777 

and the p-th iteration being: 

m 

X^n + p) = (1 - a)f(X t (n + p - 1)) + - + p - 1)) 

= / p (x*) = x\ i — 1, m 

As a result, the CML shows fixed points (x*,x*, - ■ -x*)j=i t ... tP 
of period p, which constitute one sychronized state of the same 
period. These orbits of period p constitute patterns of the CML. 

ii) To study the stability of the fixed points (x*,x*, ■ ■ • x*),, =1 ... p it is 
enough to perform it in X{ = {x\, ■ • ■ xj), because f p '(x*) has 
the same value for every fixed point x* of the p-periodic orbit. 

Let us calculate the eigenvalues of the jacobian matrix of the p-th 
iterate in that point. 

To calculate this Jacobian matrix, one must observe the following, 
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applying the chain rule: 



' dXj(n+p) \ 
s dXj(n) ) x * 



( (l-2=l a ) 



a_ 

m 



a_ 

m 



(1 - ^-a) 

V rn I 



a_ 

m 



a 
m 



a 
rn 



\ 



(1 



m— 1 
m 



f/ ( dX x (n+p-l) 
J \ x p) I <9X,(n) 



finally result: 



'aXi(n+£)\ 
> dXj(n) J x . 



™ ™ \ y rn > 



m m 



m m 



i._i.(l_ a )P \ 



m m 



V i-^(l-a) p i-i(l-a)f ••• -i + 2^(l-a)f / 



with eigenvalues 



i=i 

A = (1 - a) P Y[f'(x*) = (1 - a) p / p/ (^) multiplicity (m - 1) 

i=i 

(4) 

So, ( x*, • • • ,x*) is a stable point of the CML, and therefore, it is 
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in a stable sychronized state of period p, whenever x\ is the stable 
fixed point f p . Furthermore, as f p \x\) = f p/ (x*) j = 1, . . . , m 
all points have the same stability. 

Keep in mind that if in Theorem 1 p — 1. then the stationary sychronized 
state previously studied is recovered; and because of this, it will undergo the 
period-doubling process that will be described in what follows. 

1.3 Period doubling cascade of periodic synchronized 
states 

One would expect that if the function / p , from Theorem 1, undergoes a 
period doubling cascade, then the CML given by (j2J) shows a duplication 
cascade in the sychronized states of period p derived in Theorem 1. 

Theorem 2. The synchronized states of period p given by Theorem 1 un- 
dergo a period-doubling cascade as does f p . 

Proof 

The proof is straightforward using Theorem 1, simply using successive 
substitution of p by p ■ 2, p ■ 2 2 , . . . , p ■ 2 n , ... every time that the p- 
periodic orbit of / undergoes a period doubling bifurcation according 
to that theorem. 

Note: See figures [5] and El 
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2 r 




Figure 5: Period-4 travelling wave. CML with f(x;r) = rx(l — x), r = 
3.55464, a = 0.1 and e = 0.001 




Figure 6: Period-8 travelling wave. CML with f(x;r) = rx(l — x), r = 
3.566667, a = 0.1 and e = 0.001 
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1.4 A nonexistence theorem 

It has been proven, in Theorem 1, the existences of p-period synchronized 
states in the CML, formed by the points of p-periodic orbit of /. 
One may ask whether the n-tuple of the form 

(X 1 (n),X 2 (n),--- ,X p (n)) = {x*J{x*),...,p-\x*)) 
fP(x*)= x* j = l,---p 

that indicates that every oscillator is positioned in the successive points of 
the p-periodic orbit, generates a pattern of period p in the CML; that is 
to say, a travelling wave. Nevertheless, this presumption is false, as shown 
below. 

Theorem 3. Let {xj, x 2 , x*} be a p-periodic orbit of the function /, then 
the CML given by 

p 

Xl ( n + 1) = (1 - a)f(Xi(n + /«ra - 1)) i = 1, • • • , P 

P i=i 

does not have a p-periodic orbit of the form 

(XxM.X^n),... ,X p (n)) = (x*, f(x*), P'^x*)) 

x* being any of the p points of the p-periodic orbit. 
Proof 

The following initial conditions are taken 

(X 1 (n),X 2 (n), ■ ■ ■ ,X p (n)) = (x* l7 x* 2 , ...,x*) 
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After the first iteration, it will become 



X 1 (n + 1) = x* 2 = (1 - a)f(x{) + + ... + /(*;)) 

X 2 (n + 1) = x% = (1 - a)f(x*) + J (/(*•) + ... + /(*;)) 

^ x p (n + 1) = x; = (i - «)/(*;) + + ... + f( x ;)) 



therefore: 



x* 2 = (l-a)x* 2 + *(xt + ... + x;) 

x* = (i-«)4 + f(^ + ... + x;) 



X, 



operating it results in: 



S = f(^; + - + x;) 



from which it is deduced that 



; + + 



x l — x 2 — ... — x p 



in contradiction with x\ ^ x 2 7^ ••• 7^ x*. 

This negative result, about p-periodic waves, brings us to question the condi- 
tions under which they are produced. This study is conducted in the following 
section. 



17 



2 Analytical study of patterns in weakly cou- 
pled CML 



2.1 Travelling waves 

It is has been proven in Theorem 3 that the p-periodic orbit of the function 
f(x; r) is not inherited by the system, but it is easily observed that if a = 
then a wave of this period exists in the CML. Given that for a = the wave 
exists, one would wonder if for a small coupling a< 1, the CML admits a 
pertubative solution. For this study, we will substitute a with ea, having 
£<1 and assuming the new a is 0(1), in the CML given by (pQ). 

Theorem 4. Let {a;*, x%, x*} be a p-periodic orbit of a C 2 function /, 



such that f p '(x*) ^ 1, % 



1, . . . , p, then the CML given by 



p 



Xl ( n + 1) = (1 - ea)/(X 4 (n)) + - 



(5) 



c 3=1 
% = 1, . . . , p £ < 1 



shows a p-periodic solution given by 



Xi(n + j)=x* +j + eA. 



i+j 



% = 



..,p 



J = o, . 



p — 1 
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where 



a 



k = l,...,p 



with periodic conditions 



Ai+p — Ai 



X i+p ~ X i 



i = l, ...,p 



Proof 



The periodic orbit given by 



Xi(n + j) = x* +j + eA i+j 
i = l,...,p 
j=0,...,p-l 



will exist when the following system 



Xi(n) = x* + eAi 
Xi(n + l) = x* +1 + eA i+l 

Xi(n + p - 1) = x* +p _ 1 + eAi+p-! = x*_ x + eA^i 
Xi(n + p) = x* + eAi 



i = l,...,p (6) 



is compatible and determined. 
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As 



X^n + 1) = (1 - ea)f(Xi(n)) + - £ f( x A n )) 

^ 3=1 



from (l6l), it reults in 



* + i + sA l+1 = (1 - ea)/(xj + e^) + — ^ /(x* + eA,-) (7) 



x 



Performing the expansion 



/(x* + eAt) = /(x*) + eAif'ix*) + 0{e 2 ) 



and replacing in (JJJ), the system results in 



x* +1 + eA+i = x* i+l + eAf (xj) - eax* i+l + — ^(a$ +1 ) + 0(e 2 



i = l,...,p 



Solving the system to order e it is obtained: 



-AJ'ix*) + A i+ i = ax* +1 + 



P 



x j+i 



3=1 



i = l,...,p 



whose matricial expression is: 



o o 



-f(x* 2 ) 1 

-f'(x*) 1 



\ f Ai \ 

A 3 



2 p 3=1 J 
-^4 + ^=l X j 



... -/'(a,;) J \ A p J \ - x * + ±-X» =lX *. J 



(8) 
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This is a system of p equations and p unknowns whose coefficient matrix 
has determinant 

(-irnm) + (-i) p+1 
i=i 

Given that [T|=i f( x i) = f P '( x i) 1 (by hypothesis) the system is 
compatible and determined for every a ^ 0. Moreover, the solution of 
the system is different from the trivial one, since the independent term 
column is not null x\ ^ x* 2 ^ • • • ^ x*, that is 




It does not matter which oscillator is considered for study of the evolution 
of the system, as the algebraic system obtained is always the same. 

The solution of the system in (jHJ) can be obtained directly by inversion 
and results in: 
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a 



-f'(x{) 1 •• 

-f(x* 2 ) 1 •• 

1 0-- 








-/'(*;) ) 



X o 



1 p 



V 



/ 



The inversion of the matrix (which is not a circulant one) results in the 
following 

f Al \ 



A 2 
As 



a 



(-I)IH-I(l - (/P(xl))') 



MiV 



where the matrix M is given by 



(9) 



M : 



/'(asj) ••■/'(*;) /'(*s) ■•■/'(*;) /'W) •••/'(*;) 
i /'(xs) ■■•/'(*;)/'(*!) /'(*3) ••■/'(*;)/'(*!) 
/'(* 2 ) i /'(*;) •••/'(*;)/'(*i)/'(*5) 



i 

/'(*I)/'(*5) 



(10) 



and N by 
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1^ 



After operating in ([9j) it results in: 



A, — 2 V \ fP-j+k-i( T - 

^■k (-1)p+ 1 (1-(/p(xi))')Z^ l J 1 J+ 1 



J'=l 



k = 1, ...,j3 



Every A k ^ because the solution is known to be different from the 
trivial one. 

The solution obtained is valid at order 0(e 2 ) while s <C ^/p 7 ^*)- 



2.2 Period doubling cascade for travelling waves in a 
CML 

Period-doubling transitions to chaos have already been observed a long time 
ago, in CML with nearest neighbour coupling, using the Mandelbrot map 
[T4] . The existence of this phenomenon is not relegated only to the quadratic 
functions, and its existence can be proved for any function (as we will demon- 
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strate) undergoing a period-doubling cascade; therefore, this phenomenon 
must be very frequent. 

Theorem 5. Let /:/—>/ be a C 2 funtion depending on some parameter, 
in function of which the 2 p -periodic orbit of the map x n+ \ = f(x n ) 
undergoes a period-doubling cascade. Let |^* 2 p+9j. be the 2 p+q - 
period orbit of the cascade, q G N, where it is noted that f k (x* 2P+q ) = 

X i+k,2P+i- 

The CML given by 



2 p 

X t (n +1) = (1 - + f 

i=i 

i = l,-- ,2? e < 1 
has a 2 p+9 -periodic solution given by 

Xi(n + j) = ^29(j-l)+l+i,2P+"J + £ ^29(i-l)+l+j 

i= l,...,2f j=0,...,2^ 

where 



( _ 1+(/ 2P+« (a . 

j'=i 



2P+1 

— \j 2P q 3+k l {x* +12P+q )^ [(x 
k = 1,...,2p+« 



(12) 
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with 



' 2 p 

i=i 
2 p 

E 

i=i 

E 

i=i 

E 

i=i 



2 < J(i-l)+2,2P+9 



29(i-l)+3,2P+9 



x; 



2'3(i-l)+4,2P+9 



2<J(i-l)+5,2P+9 



2 p 

E 

i=i 



:r 



29(i-l)+2«+l,2P+« 



if j = [2 9 ] 

if j = [2«] + 1 

if j = [2 q ] + 2 

if j = [2 9 ] + 3 

if j = [2 q ] +2 9 -l 



where [2 q ] represents a multiple of 2 q . This 2 p+l3 -periodic solution fulfills 
the periodicity condition 



Proof 



Note: Notice that although there are just 2 P oscillators the periodic 
wave will have period 2 p+q , therefore, the system that is dealt with in 
the proof will have 2 p+q equations. 

Let 2 P oscillators be with the initial conditions given by 



Xi(n) = x* 2q(i _ 1)+h2P+q + eA 2q{i _ l)+l i = 1, • • • , 2 P (13) 



Initial conditions are fixed points of f 2P+9 (taken one every 2 q ) plus a 
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perturbation that must be calculated. A 2 p+9 -periodic orbit will exist 
whenever the system 

Xi(n) = a ; 2«(i-i)+i,2P+i + e ^2<j(i-i)+i 

Xi(n + 1) = x 2i(i-l)+2,2 p+q e ^2«(i-l)+2 

i i = l,...,2" (14) 

Xi{n + 2 p+q ) = 2 ; 29(i-l) + l+2P+i,2P+9 £ ^ L 2i(i~l)+2P+i 
= X 29{i-l)+l,2"+i EV ^29(i-l)+l = ^i( n ) 

is compatible and determined. 

From equation f fT2l) and subtituting the second equality in f [T4l) . we 
have 

X 29(i-l)+2,2P+<2 + £ ^2?(i-l)+2 = (1 — ^Oi) f { x *2i(;i-l)+l,2P+i + £ ^29(i-l)+l) + 

ae 

+ 2P Z^^*- X 29(j-l) + l,2P+9 + £A 24(j-l)+l) 
3=1 

(15) 

Performing a Taylor expansion of / to order 0(e 2 ) and substituting in 
(1151) the following is obtained 

X 2<i(i~l)+2,2P+i e ^2i(i-l)+2 = X 2i(i-l)+2,2P+i^~ 
+£y4.2?(i— 1)4-1 / ( a '2 | 3(i-l)+l,2P+'j) ~~ £Q; * E 29(j-l)+2,2P+<2 2? X/i=l 2 '2'J(i-l)+2,2J > +9 ) 

Doing exactly the same with the next equality in <HM the following is 
obtained: 

X 29(j-l)+3,2P+9 + £ ^29(j-l)+3 — X 2i(i-l)+3,2P+i~\~ 

2P 

+£y4 2 9(i-l)+2/'(2 ; 2'J(j-l)+2,2P+'j) — £aX 2i(i-l)+3,2P+i + T^T 5Z a '2<?(i-l)+3,2P+9 + ^K^) 

i=l 
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and with the last equality in (fill) , we get the following equation: 



2«(i-l)+ 



1 2P+1 £ ^-2«(i-l)+l 



X. 



+eA 2 < 1 (i-i)f'(x* 2q(i _ 1 



29(i-l)+l,2P+« 

ea 



+ 



(i-X),2P+i) £aX 2i(i-l)+l 



i=l 



The former 2 p+(? -equations, for the oscillator i, represent a linear system, 
whose matricial expresion is: 



/ (^(i-ij+i^p+g) 



1 ••• 



* / ^2«(i-l)+l ^ 
^29(i-l)+2 



l)+2,2P+9 ' 2P 2^ X 2l(i-1)+2,2P+1 



/'( X 29(i-l) + 2P+ 1 ,2P+9) / V ^29(i-l) + 2P+9 / 
2 P+1 



2l(i- 



i=l 
2 P+ 



29(i-l)+3,2P+9 T 2P 



)+3,2P+9 



29(i-l) + l,2P+9 2 P 
\ i — 1 



op+1 

i v » 



2p+g 



(16) 



The previous matricial expresion represents a system of 2 p+g equations 
with 2 p+q unknowns, and being the determinant of coefficient matrix 



2P+1 



-i? p+q n /'< 



x*) - (-1 



2P+1 



8=1 



-MO 



(period doubling bifurcations take place when f 2P+9 {x*) = — 1). 

Thus the system is compatible and determined for every a and, as in 
the previous cases, its solution is different from the trivial one for a^O. 
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The solution is obtained directly from (PT6l) by inversion: 



( A \ 



.4 



29(i-l)+2 



y ^42<?(j-i)+2P+9 J 



(-l+[f 2P+q {xl)]') 



MN 



(17) 



where 



M 



1 -n 





••• 

1 








-1 



X 2i(i~l)+2P+i,2P+i) ) 



has already been calculated (see [TQJ) and 



/ 



N 



2PH 



\ 



29(i-l)+2,2P+9 n 2P / ■> 2gCi-ll+2.2P+g 

i=l 
2 p+i 



28(i-l)+3,22>+8 ~ 2P / / 2g(i-l)+3,2P+« 
i=l 



2/H 



■' •->„ , i , i ■■>,. i-o + _ 



2«(i-l)+l,2P+« ^ 2P / / 29^-11+1. 2P+ 1 
i=l 



After operating in (11711 it results in: 



Ai 



a 



(-1)^(1 -[p p+q (xW)u 



2P+1 

^2 \ f w+*- j+ k-i (x . 

r- 



j+l,2P+i) ( { X l+j,2P+i) + ^3 
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with 



k = 1, 2 p+<? 



' 2 p 

i=i 
2 p 

E 

i=i 

E 

i=i 

E 



29(i-l)+2,2P+9 



2'3(i-l)+3,2P+9 



29(i-l)+4,2P+9 



29(i-l)+5,2P+9 



2 p 

E 



x 



28(i-l)+2«+l,2P+9 



if j = [2«] 

if j = [2«] + 1 

if j = [2 q ] + 2 

if j = [2 9 ] + 3 

if j = [2 q ] +2 9 -l 



where [2 q ] represents a multiple of 2 q , with the periodicity conditions 



i+2P+i 



= As 



-V' — — / ' 



due to the cyclic character of the 2 p+9 -periodic orbit. 

The solution obtained is valid at order 0(e 2 ) while e <C — p + g — . 

i — [/ ( x i)\ 

The saddle-node orbit ([/ 2P+ "(^)] / ^ 1) has been avoided by the con- 
dition q G N (0 ^ N), and therefore there has been, at least, one 
period-doubling bifurcation. 
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Remarks: 

1. Notice that this theorem indicates that a CML, with 2 P oscillators, has 
originally a 2 p -period travelling wave. As the 2 p -periodic orbit of / 
duplicates (q increases) to a 2 p+(? -periodic orbit so does the travelling 
wave of the CML. 

2. The theorem 5 does not impose any restriction to the 2 p -periodic orbit 
of /. In the case of / presenting 2 p -periodic windows, this 2 p -periodic 
orbit could belong to a period-doubling cascade in the canonical win- 
dow, or originate from a 2 p -periodic saddle-node orbit. In the former 
case, / undergoes a period-doubling cascade, in the latter it is f 2P 
who goes through a period-doubling cascade (this period-doubling cas- 
cade would be located inside a 2 p -periodic window). The conclusion 
is straightforward: the CML will not have just one 2 p -periodic wave, 
undergoing a period-doubling cascade; there will be as many as 2 Pl - 
periodic windows of /, with p 1 < p (p 1 = would be the canonical 
window) . 

3. Similar arguments can be done for synchronized state cascades, subject 
to the condition that the duplicating orbit does not have prime period: 
it can be a p ■ (/-periodic orbit in the canonical window, or a (/-periodic 
orbit in the p-periodic window, or a p-periodic orbit in a (/-periodic 
window, that afterwards will undergo period bifurcation cascade. 

4. There is a fact that could be not observed at a first reading of theorem 
5: the points used to construct the perturbative solution both can be 
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stable and unstable. 

5. Since it has been deduced that the CML undergoes a period doubling 
cascade and that this cascade has its origin in the period doubling 
cascade of the /, it is concluded that the CML inherits the dynamics 
off. 

3 Discussion and conclusions 

Several theorems have been proved that show the sufficiency conditions of 
existence of synchronized states (periodic and chaotic) and travelling waves 
in CML. Also it has been analytically determined the value that describes 
the state of each oscillator at any moment. The results of the theorems are 
as general as possible. This is due to two facts. Firstly, the CML, with 
which we have worked, has a number of arbitrary oscillators. Second, the 
function /, that rules the dynamics of every oscillator, is also arbitrary with 
the condition that it undergoes a period-doubling cascade. 

The results have the following consequences and link with other research: 

i) The emergence of the global properties from the local ones has been 
proved. The global dynamics inherits the dynamics of every oscillator: 
fixed points of the system come essencially from the fixed points of the 
map (that governs the dynamics of every oscillator) compounded with 
itself m-2 k times (being m number of oscillators of the CML). In partic- 
ular, this result has been observed recently in numerical computations 
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Our results are an explicit analytical expresion of the results of Lemaitre 
y Chate [TH], who proved, in CML, the traslation of the local properties 
to a spatiotemporal level. 

) The dynamics of a CML has been studied, where the individual dy- 
namics of every oscillator is ruled by an arbitrary function /; being 
few the analytical results on the matter, one normally only works with 
quadratic functions or piece- wise linear functions [TUl [T2l E] ■ The one 
presented here is an interesting generalization that permits the calcu- 
lation of properties associated with the states and their evolution. 

) Two limitations that are present when numerical techniques are applied 
have been overcome: 

(a) Spurious results, due to finite precision in computer simulations 
(HUH [20], have been avoided. 

(b) Limits tending to infinity can be used for analytical solutions, 
both with the number of oscillators in CML and the number of 
bifurcations in the period doubling cascade. 

On the one hand, the numerical simulation with a large number 
of oscillators becomes unaffordable due to the computation time 
that would be necessary as the number of oscillators grows. On the 
other hand, as it has been indicated in the introduction, the study 
of the onset of the turbulence in a fluid to be properly understood 
would need many oscillators, the more the better. 

As a direct consequence, of taking the limit in the period doubling 
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cascade, it is deduced the existence of waves of arbitrary period, 
tending to infinity as the parameter bifurcation gets closer to the 
Myrberg-Feigenbaum point. This is a response to the established 
question of Gade and Amritkar in their work [21] where they found 
the wavelength-doubling bifurcation. Another question raised by 
Gade and Amritkar in that same paper was: is there more than 
one value such that, if the parameter value tends to it, then the pe- 
riod of the travelling wave tends to infinity?. The response again 
in the afirmative; in fact, there are infinite values that are the 
correspondent Myrberg-Feigenbaum points of the windows inside 
the canonical window. The position of these values is determined 
by the Saddle-Node Bifurcation Cascades [22] and the relation be- 
tween period doubling cascade and Saddle-Node Bifurcation Cas- 
cades is also known [23j. 
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